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Abstract 



Pl,: 

6 . 

£3 , We determine numerically the parallel, perpendicular, and antisymmet- 

ric diffusion coefficients for charged particles propagating in highly turbulent 
magnetic fields, by means of extensive Monte Carlo simulations. We propose 
simple expressions, given in terms of a small set of fitting parameters, to 
account for the diffusion coefficients as functions of magnetic rigidity and tur- 
bulence level, and corresponding to different kinds of turbulence spectra. The 
results obtained satisfy scaling relations, which make them useful for describ- 
ing the cosmic ray origin and transport in a variety of different astrophysical 
environments. 

The diffusion and drift of charged particles across highly turbulent magnetic fields 
are key issues in describing the transport of cosmic rays in different astrophysical 
environments, e.g. the interplanetary, interstellar and intergalactic media, as well 
as the efficiency of Fermi acceleration processes at cosmic ray sources. In particular, 
it has been shown that the inclusion of drift effects in the transport equation leads 
naturally to an explanation for the knee, for the second knee and for the observed 
behavior of the composition and anisotropies between the knee and the ankle [1, 2, 
3, 4]. However, the accuracy of the investigations performed so far are limited by 
the lack of conclusive results concerning the behavior of the diffusion tensor under 
highly turbulent conditions as a function of the particle energy and the relevant 
magnetic field parameters. In particular, the magnetic fields in the Galaxy are 
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highly turbulent because the mean random field is of the order of the mean regular 
field. Moreover, since there are reversals in the orientation of the regular field, this 
implies the existence of regions with negligible regular fields in which the turbulence 
prevails. 

Perturbative studies for low turbulence have been developed since long ago [5, 
6, 7], but these analytic methods cease to be applicable for high turbulence levels, 
and only recently the parallel and transverse diffusion coefficients were calculated 
numerically for regimes with high turbulence [8, 9]. The aim of this work is to 
provide a thorough and more systematic calculation of these coefficients, and to 
parametrize the results in order to make them useful in a variety of different kinds 
of applications. Moreover, we present here also a numerical evaluation of the Hall 
diffusion coefficient that is responsible for the drift effects, which so far has never 
been evaluated quantitatively under highly turbulent conditions. It should also be 
remarked that, while in [8, 9] only the Kolmogorov spectrum of fluctuations in the 
random magnetic field was considered, in this work other types of turbulence are 
studied as well (namely, the Kraichnan and Bykov-Toptygin turbulence spectra, 
which bracket a wide range of possible turbulence spectra). 

A relativistic particle of charge Ze propagating in an uniform regular magnetic 
field B describes a helical path characterized by a pitch angle 6 and a Larmor 
radius given by 

" ZeB ~ 10 15 eV \/iG ) P ' 1 ] 

The component of the velocity parallel to B is v» = ccos#, while the radius of 
the helical trajectory is r^sin^. In the presence of a random magnetic field B r 
with a maximum scale of turbulence L max , the particles scatter off the magnetic 
irregularities and change their pitch angle, but not their velocity. The pitch angle 
scattering proceeds mainly in resonance (i.e., the scattering is dominated by the 
inhomogeneities with scales of the order of r£), and hence it is an effective mechanism 
of isotropization as long as < L max . For instance, for the galactic magnetic field, 
with strength B ~ few /iG and maximum scale of turbulence L max ~ 100 pc, the 
pitch angle scattering leads to a diffusive regime for protons with energies up to 
few 10 17 eV. 

In general, the diffusion tensor can be written as 

Dij = (D\\ - D ± ) M, + DJij + D A e ijk b k (2) 

where b = B /-Bo is a unit vector along the regular magnetic field, 5ij is the Kro- 
necker delta symbol, and is the Levi-Civita fully antisymmetric tensor. The 
symmetric terms contain the diffusion coefficients parallel and perpendicular to the 
mean field, D\\ and D±, which describe diffusion due to small-scale turbulence, while 
the antisymmetric term contains the Hall diffusion coefficient D^. 
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The diffusion along the magnetic field direction is due to the pitch angle scat- 
tering and leads to a diffusion coefficient given by 



D \\ = f A,, , (3) 

where Ay is the mean free path in the parallel direction [1]. In this expression, Ay 
depends on the power of the random magnetic field modes at scales of order of re, 
i.e. 



An oc 



dE r /d Ink 



(4) 



where dE r /dk is the power spectrum of the random magnetic field energy density. 

In general, a component of the diffusion tensor is associated to some particular 
physical processes which may not contribute to the other components. The diffusion 
transverse to the regular magnetic field is due to both pitch angle scattering, which 
is the mechanism that prevails in the parallel diffusion, and to the wandering of 
the magnetic field lines themselves, which drag with them the diffusing particles in 
the direction perpendicular to Bo [6, 10, 11]. When the turbulence level is small, 
the diffusion in the orthogonal direction is much slower than in the parallel one 
and is strongly affected by the field line random walk, but in the limit of very high 
turbulence, the parallel and perpendicular motions become similar. 

The Hall coefficient Da describes in turn the macroscopic drift associated to the 
gradient of the CR density, leading to a current 

J A = D A b x X7N . (5) 

This macroscopic drift is orthogonal to both Bo and \7N, and in particular is also 
present for a constant regular field. The relation between this macroscopic current 
and the microscopic drift associated to gradients and curvature in B has led to some 
confusion and controversy in the literature [12, 13, 14, 15]. A common approach is 
to identify a guiding center, associated to the instantaneous radius of curvature 
of the particle's trajectory, and to define a guiding center velocity v g . In the so- 
called first order orbit theory, an average over a particle's gyroperiod is performed 
assuming that the scale of variations in the field is much smaller than the Larmor 
radius, and then an ensemble average over a given distribution of particles is carried 
out. Alternatively, (v g ) can be calculated by averaging locally over a distribution of 
particles in phase space [15]. In the latter case, the results are given in terms of an 
expansion of the anisotropy of the momenta in spherical harmonics, but in principle 
no assumptions concerning the scale of variations of the field are required. On the 
other hand, the average particle velocity (v) can be calculated from the Vlasov (or 
collisionless Boltzmann) equation, averaging again over a particle distribution in 
phase space. As shown in [15], the mean particle and guiding center velocities are 
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related through their definition, and they differ in a term which is sometimes called 
diamagnetic drift. A qualitative difference that can be pointed out is that (v) has a 
contribution arising from density gradients and can be nonvanishing even in the case 
of an uniform regular magnetic field, while (v g ) depends only on the field gradients. 

Notice that, when making an average of the particle velocities inside a given 
volume element, one is including particles whose center of gyration are outside the 
volume considered, and, inversely, considering the guiding centers inside a given 
volume may correspond to taking into account particles which are actually outside 
that volume element. Since the observable quantity is the average particle velocity 
(v) and not the abstract concept of the guiding center motion (v g ), it is clear that 
it is the former one that will enter into the drift motions. 

As commented above, the diffusion coefficients were calculated analytically only 
in the case of small turbulence levels, while results valid under highly turbulent 
conditions require instead a numerical approach [8, 9]. Following [8], we will con- 
sider charged particles propagating in a magnetic field of the form B(r) = B z + 
B r (r), where the first term represents an uniform regular field directed along the 
z— direction, while the second corresponds to the random component. In order to 
approximate numerically the isotropic and spatially homogeneous turbulent field, 
one can sum over a large number (N m ) of plane waves with wave vector direction, 
polarization and phase chosen randomly [16, 8], i.e. 

N m 2 

BrW = E E Mk n )C cos(k„ • r + C) , (6) 

n=l a=l 

where the two orthogonal polarizations £° (a = 1,2) are in the plane perpendic- 
ular to the wave vector direction (i.e., _L k n , so as to ensure that V • B = 0). 
The wavenumber distribution is taken according to a constant logarithmic spac- 
ing between k min = 2ir/L max and k max = 2n/L min , where L min and L max are the 
minimum and maximum scales of turbulence, respectively. The energy density of 
the turbulent component is taken as dE r /dk oc A; -7 , where the spectral index 7 is 
given by the kind of mechanism that builds up the turbulence. Hence, the plane 
wave amplitudes satisfy A 2 (k n ) = J\f(B 2 )k~' r (k n — fc n _i), with M a normalization 
constant that assures that J2n A 2 {k n ) = (B 2 ). In this work we will consider in par- 
ticular three spectra of interest for astrophysical applications, namely a Kolmogorov 
spectrum with 7 = 5/3 (which is the only case studied numerically in the past, and 
is particularly attractive since, according to observations [17, 18], the density fluc- 
tuations in the interstellar medium follow this turbulence spectrum), a Kraichnan 
hydromagnetic spectrum with 7 = 3/2 [19], and the Bykov-Toptygin spectrum with 
7 = 2 [20]. 

By means of the Kubo formalism [21, 22, 23], the diffusion coefficients can be 
computed directly by taking ensemble averages of the decorrelation between different 
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components of the single particle velocities, i.e. 



oc 



Dij = j dtRij(t), (7) 

with Rij(t) = (voiVj(t)) (with i> j = Vi(t = 0) and where (...) denotes the ensemble 
average taken over an isotropic distribution of many particles). In [23], it was 
assumed that the velocity decorrelations were modulated by exponential factors, 
adopting the following ansatze 

c 2 

Rxx(t) = Ryyit) = - cos cot e~ t/T± , (8) 
o 

c 2 

R yx (t) = -R xy (t) = - sin cut e' t/TA (9) 

and 

R zz (t) = j e-«/n , (10) 

where u> = c/r L is the Larmor angular gyrofrequency, and where t\\,t± and ta are 
the decorrelation timescales associated to the different diffusion components. Then, 
the diffusion coefficients were obtained by integrating the proposed expressions for 
Rijit) using Eq.(7). The diffusion coefficients were found to be 

D\\ = , (11) 

D ± = ^ ^ 2 (12) 
3 l + (^r ± ) 2 1 J 



Da = ^t^—, ■ (13) 



and 

cr L {u)T A y 

3 \ + {ujT A y 

It is interesting to note that expressions of this form, but in which only a single 
timescale r appears for the three diffusion coefficients, were obtained also in other 
analytic approaches that assumed a single scattering process to be responsible for 
all the decorrelations [24, 25, 26]. However, it should be pointed out that the 
expressions proposed in [23] for Rijit) assume implicitly a small departure from the 
helical trajectories, and they are no longer adequate for high turbulence levels (see 
below for further discussions). Moreover, there is no general theory providing the 
decorrelation timescales, and this requires then further assumptions in order to make 
Eqs.(ll)-(13) useful. These additional assumptions are sometimes obscure and lead 
to results often at odds with the outcome of numerical simulations, as pointed out 
in [2]. 
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Alternatively, the parallel and perpendicular diffusion coefficients can be calcu- 
lated from the asymptotic rate of increase of the mean squared displacements in 
each direction, namely 

D , = lim «M!> (14) 
11 At^oo 2At 

D L = lim ^fl = Urn M . (15) 

At^oo 2At At^oo 2 At 

In any case, the straightforward method consists basically in generating a sample 
configuration for the random magnetic field (by choosing randomly the propagation 
direction, polarization and phase of the N m plane waves in Eq.(6)), and then follow- 
ing the trajectory of a particle that propagates in the total (regular plus random) 
field from the origin with a random initial direction. The results should then be av- 
eraged over a large number of different field configurations in order to calculate the 
corresponding diffusion coefficients. In [8], for instance, the results were obtained 
by integrating firstly the trajectories of 2500 particles for a given field configuration, 
and then by averaging over 50 different field realizations. We found, however, that 
the results are more biased by the particular choice of field configuration than by 
the choice of the initial velocity of the particle. Hence, it turns out to be more 
convenient to average directly over a large number of field realizations. Indeed, we 
typically generated ~ 10 5 field configurations, following the trajectory of a single 
particle in each of them, and used N m = 100 modes in all simulations. 

For the numerical integration of the particle trajectories under the influence of 
the Lorentz force, we adopted a time step At = 0.1 r L /c in the Runge-Kutta rou- 
tine. In [9], the parallel and perpendicular diffusion coefficients were computed 
from Eqs.(14)-(15), and this method requires to follow the particle trajectories 
for quite long times in order to reach the asymptotic region, typically longer than 
t = 1000 tl/c. We found more convenient instead to compute directly the particle 
decorrelation functions Rij(t) and integrate them through Eq.(7), since this proce- 
dure requires to follow the particle trajectories for times typically not larger than 
t = 100 tl/c. The underlying reason for this seems to be that the displacements 
keep more memory of the initial velocity adopted than the velocities themselves. In 
addition, this method allows to compute the antisymmetric diffusion coefficient by 
means of R yx , while averages such as {Ax Ay) vanish for large times, and hence they 
are not useful for computing Da- 

Figure 1 shows the time dependence of the decorrelation functions associated 
to different velocity components, for the Kolmogorov case and different turbulence 
levels, with the turbulence level defined here 1 by a 2 = (B 2 )/Bq. For low turbulence, 

1 Alternatively, one can define the turbulence level as 77 = (B 2 )/({B 2 ) + B 2 ) = cr 2 /(l + a 2 ) [9]. 
The definition adopted here follows that in [8]. 
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Figure 1: Time dependence of the particle decorrelation functions associated to 
different velocity components, for the Kolmogorov turbulence spectrum and two 
different turbulence levels: (a) a 2 = 0.3 and (b) a 2 = 5, both corresponding to 
tl — 0.1 L max . 
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Figure 2: Linear-log plots of the local maxima of \R XX \ and \R yx \ versus t corre- 
sponding to the case of Fig. 1(a), showing that tj_ = ta (see Eqs.(8)-(9)). The lines 
are guides to the eye. 



it is observed in Figure 1(a) that, after a given initial transient period, the correla- 
tion functions behave as Eqs.(8)-(10), and hence one can calculate the decorrelation 
timescales directly. For instance, following the local maxima of the sinusoidal func- 
tions one finds t± = r A , as it is apparent from Figure 2. In [23], these decorrelation 
timescales were considered to be equal as a simplifying assumption, but it was sug- 
gested that new effects could arise in a general case with r± ^ ta- Here we find that 
these decorrelation timescales are indeed equal, and actually this could be expected 
from the fact that the physical origin for the decorrelations is the same in both cases. 
For high turbulence levels, however, the amplitude decrease is much more abrupt, 
as it is shown in Figure 1(b). The decorrelation functions become vanishingly small 
already in the transient period, and Eqs.(8)-(10) are hence not valid anymore. 

Defining the dimensionless parameter p = r^j 'L max , the results given in terms of 
Dj (cL max p) vs. p are universal and can be scaled in order to calculate the diffusion 
tensor for different sets of values for the regular field amplitude, random field length 
scale and particle energies. For instance, the range 0.01 < p < 1 investigated in 
this work can be regarded as corresponding to 10 15 < E/eV < 10 17 for protons 
propagating in the Galaxy (for Bq = 1 pG and L max = 100 pc). Alternatively, 
considering protons propagating in the interplanetary field, with strength Bq = 



8 



0.01 



0.1 1 

p 



10 



Figure 3: Results obtained for the parallel diffusion and corresponding to the Kraich- 
nan spectrum for different turbulence levels in the range O.f < a 2 < fOO, as indi- 
cated. The lines are guides to the eye. 

50 pG and L max = 0.01 AU, the limit to the diffusion regime at p = 1 corresponds 
to the kinetic energy Ef~ = 1.8 GeV. 2 Since pitch angle scattering proceeds mainly in 
resonance, the particle propagation is essentially independent of the minimum scale 
length of turbulence L min , as long as L min C r L . For definiteness, we then adopt 

L m in = 0.1 Tl- 

Figure 3 shows the numerical results obtained for the parallel diffusion coefficient 
corresponding to the Kraichnan spectrum and for different turbulence levels in the 
range 0.1 < o 2 < 100, as indicated. As follows from Eqs.(3)-(4), in the resonant 
scattering regime one expects that D\\ oc p 2 ~ 7 , while at high rigidities (outside the 
resonant scattering regime) the scattering effectiveness is expected to decrease as 
E~ 2 [1, 27], hence leading to D\\ oc p 2 , irrespective of the turbulence spectrum con- 
sidered. The behavior of the results displayed in Figure 3 agrees very well with these 
expectations, showing a crossover between the resonant and non-resonant scatter- 
ing regimes that takes place around p ~ 0.2. Indeed, the actual scale length that 
separates both regimes is determined by the correlation length of the turbulence 

2 In [8], however, the power spectrum of the random interplanetary field is considered to con- 
tinue with constant amplitude from 0.01 AU up to a maximum scale of 1 AU, hence extending 
considerably the diffusive regime to higher energies. 
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Table 1: Parameters of the fitting formulae corresponding to the parallel, trans- 
verse and antisymmetric diffusion coefficients, and for different kinds of turbulence 
spectra. See more details in the text. 



spectrum L c , defined as 

/oo 
dL(B r (0) • B r (r(L))> = L c (B 2 r ) , (16) 
-oo 

where the point r(L) is displaced with respect to the origin by a distance L along a 
fixed direction. Then, considering a spectrum of fluctuations of index 7 extending 
between the scale lengths L min and L max , the correlation length is given by [28] 

j 1 t 7 1 J {J-'min I Lmax) ( 1 7\ 

~ o, ax ^ 1 _ f /" 77 vt^T • v 1 ' J 

7-1 \J->min I J-'max ) 

Hence, L c /L max ~ 0.2 is a quite representative value for the random field power 
spectra considered in this work, and this explains the change of regime observed at 
p ~ 0.2. 

According to these considerations, an appropriate way of fitting the results is to 
interpolate Dn between the power laws that characterize the low- and high-rigidity 
regimes. A convenient way to achieve this is by means of the expression 



D\\ N\\ 
cL rnax p a 2 \ 




(18) 



where the parameters N» and p» are given in Table 1. This expression appears to 
fit well our numerical results up to a 2 ~ 10, which is already a very high turbulence 
level in most astrophysical applications of interest. For even higher turbulence levels, 
the results are best described asymptotically by the parameters N^ smt and p|| Smi , 
also given in Table 1. 

In order to check the consistency with the previous numerical results for highly 
turbulent parallel diffusion (which was studied in [9] assuming a Kolmogorov spec- 
trum of fluctuations), Figure 4 shows a comparison between our fit and the re- 
sults of [9], both corresponding to a Kolmogorov spectrum and for the same tur- 
bulence levels. Notice that in [9] the Larmor radius 3 r* L is defined by replacing 

3 To avoid any confusion, we adopt an asterisk to refer to the quantities denned in [9]. 
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Figure 4: Comparison between the fit to Dn given in this work (formally only valid in 
the range 0.01 < p < 1) and the results of [9], both corresponding to a Kolmogorov 
spectrum and for the same turbulence levels. 



Bq + B 2 in Eq.(l), hence coupling in r* L the dependence on rigidity and on 



y "u 1 "r 

the turbulence level 



while, on the other hand, the dimensionless parameter p* is 
Hence, the results have to be rescaled according to the 



defined as p* = 2nr* L / L m 
relations D/(cL max p) = D/{r* L c)/y/\ + a 2 and p = p*\J\ + (t 2 /2tt. As can be seen 
in Figure 4, both data sets agree reasonably well. 

The perpendicular diffusion coefficient is most easily parametrized from the re- 
sults for the ratio D±/D», since this ratio exhibits little dependence with the rigidity. 
In order to avoid the subdiffusive regime (see discussion below), we first restricted 
ourselves to p > 0.03 and a 2 > 1. As in previous numerical investigations [8, 9], 
we found the ratio D±/D\\ to be independent of rigidity in the low-rigidity region 
(namely, for p < 0.2), while it scales as p~ 2 in the high-rigidity region. Hence, our 
results can be conveniently parametrized by means of the expression 



x(^) 



2\a ± 



X 



1 (p < 0.2) 

(p/0.2r 2 (p > 0.2) 



(19) 



where the parameters N± and a± are also given in Table 1. In [2, 3, 4], a fit to the 
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Figure 5: Fit to D±/D\\ as a function of rigidity for the Kolmogorov spectrum 
and corresponding to the low-rigidity regime (i.e. p < 0.2), as obtained for different 
turbulence levels. For larger rigidities, D±/D\\ oc p~ 2 (see Eq.(19)). For comparison, 
results obtained from fits to the data given in [8] and [9] are also shown. 



data of [9] (for the low-rigidity regime) was given as 



D 



_L 



1 



D\\ 1 + 4.5 2 /( ( t 2 ) L5 ' 



(20) 



while the results presented in [8], which correspond to the turbulence range 0.03 < 
a 2 < 1, can be accounted for by the expression 



= 0.025 x (a 2 ) 1 ' 835 . 



[21] 



Figure 5 shows the fit given by Eq.(19) for the Kolmogorov case and different 
turbulence levels. For the sake of comparison, the results of [9] (as given by Eq.(20)) 
and of [8] (as given by Eq.(21)) are also shown. For moderate turbulence levels 
(a 2 = 1) our results show an excellent agreement with the results of [8], while the 
agreement with [9] tends to be better for larger turbulence. 

For low rigidities and turbulence levels (for instance, p ~ 0.01 and a 2 < 1), 
we have found evidence for the phenomenon of subdiffusion, as already reported in 
[9]. Indeed, a low-rigidity particle at low turbulence tends to remain attached to a 
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Figure 6: ((Ax^) 2 ) /2At/ \cL rnax p) versus time for different values of rigidity and 
turbulence levels, in the case of a Kolmogorov spectrum of fluctuations. The subd- 
iffusive regime shows up at low enough rigidities and turbulence levels. 

field line, and hence transverse diffusion chiefly proceeds by the transverse random 
walk of field lines. The mean perpendicular displacement is then expected to evolve 
more slowly with time than in the case of normal diffusion (i.e., (Axj_) oc t m with 
m < 1). Figure 6 shows ((Ax ±) 2 ) /2At / (cL max p) as a function of time, for different 
values of rigidity and turbulence levels, in the case of a Kolmogorov spectrum of 
fluctuations. While a plateau is attained at large t for p = 0.1 and a 2 = 1, hence 
corresponding to the usual diffusion relation (Ax 2 ± ) oc t, for lower values of rigidity 
and turbulence level the phenomenon of subdiffusion shows up. As an approach to 
understanding the subdiffusive regime, it has been proposed and investigated the 
so-called compound diffusion, in which particles are assumed to be strictly tied to 
the field lines, while they scatter back and forth along the lines [29, 22, 30, 31]. For 
the limiting case of compound diffusion, it turns out that m = 1/2, while in the case 
of three-dimensional particle transport m is expected to have a smooth dependence 
with rigidity and turbulence, such that 1/2 < m(p,a 2 ) < 1. Further investigations 
aiming at a detailed, quantitative description of subdiffusion are currently under 
progress. 

Concerning the antisymmetric diffusion coefficient, the numerical results can be 
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Figure 7: Comparison between the fit to D A given in this work (for the Kolmogorov 
spectrum and the rigidities p = 0.01, 0.1, 1) and the simple p— independent prescrip- 
tion adopted in [2, 3, 4] (see Eq.(24)). 



fitted by the expression 

Da _ 1 1 

cL max p 3 + (a 2 /alf 

where 



(22) 



^=^ X {L9p-(P>0.2 2 j ' < 23 » 

and where the values for the parameter N A are given in Table I. In the limit of 
very low turbulence, this expression tends to the appropriate value Da — cr L /3 [1] 
(see also Eq.(13) in the large Ta limit), while it vanishes in the limit of very strong 
turbulence, as expected. In [2, 3, 4], the propagation of cosmic rays diffusing in the 
Galaxy was studied and, due to the lack of an expression like that given by Eqs.(22)- 
(23) (i.e. valid even under highly turbulent conditions), a simple p— independent 
prescription for Da/ cL max p was adopted, namely 

DA 1 1 (24) 



cL max p 3 l + (a 2 ) L5 /4.5 2 ' 
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which was inferred for a Kolmogorov spectrum of fluctuations. In Figure 7 this 
ad hoc prescription is compared to the actual fit of Eqs.(22)-(23), calculated for 
the Kolmogorov case and the rigidities p = 0.01,0.1, 1. It can be observed a sound 
agreement between the prescription of [2, 3, 4] and the fit corresponding to p = 0.01, 
and, in any case, the dependence of D^/p with p is apparent only for very high 
turbulence levels a 2 ^> 1. 

As a summary, in this work we performed extensive Monte Carlo simulations 
to determine the parallel, transverse and antisymmetric diffusion coefficients that 
describe the propagation of cosmic rays under highly turbulent conditions. We 
examined the simple analytical approach proposed in [23], and found that the ex- 
pressions given there in terms of mean trajectory decorrelations are only meaningful 
for low turbulence levels. Furthermore, we found that, as long as their approach is 
valid, the decorrelation timescales associated to transverse and antisymmetric dif- 
fusion are indeed equal. We evaluated the diffusion coefficients and parametrized 
the results by means of simple expressions, which agree with the expected behavior 
in the limit of low turbulence levels. Moreover, the results obtained were com- 
pared to the previous numerical calculations performed in [8, 9]. In this respect, 
this work extends the previous investigations, since it also takes into consideration 
other possible turbulence spectra in addition to the Kolmogorov case (namely, the 
Kraichnan and the Bykov-Toptygin spectra), it provides useful parametrization for- 
mulae, and it includes the study of the antisymmetric diffusion coefficient. Finally, 
we compared the new results for the antisymmetric coefficient with the prescription 
adopted previously in [2, 3, 4] for explaining the cosmic ray spectrum, composition 
and anisotropies in the region between the knee and the ankle. We hope that these 
results will be useful in a variety of different astrophysical scenarios related to the 
origin and transport of cosmic rays. 

Acknowledgments 

Work partially supported by CONICET and Fundacion Antorchas, Argentina. E.R. 
is partially supported by a John Simon Guggenheim Foundation fellowship. J.C. is 
currently supported by the Program for Latin American Students of the Theoretical 
Physics Department of Fermilab. Fermilab is operated by Universities Research 
Association Inc. under contract no. DE-AC02-76CH02000 with the DOE. 

References 

[1] V.S. Ptuskin et al, Astron. Astrophys. 268 (1993) 726. 

[2] J. Candia, E. Roulet and L.N. Epele, JHEP 0212 (2002) 033. 



15 



[3] J. Candia, S. Mollerach and E. Roulet, JHEP 0212 (2002) 032. 
[4] J. Candia, S. Mollerach and E. Roulet, JCAP 0305 (2003) 003. 
[5] E.N. Parker, Planet. Space Sci. 13 (1965) 9. 
[6] J.R. Jokipii, ApJ 146 (1966) 480. 

[7] S. Chapman and T.G. Cowling, The mathematical theory of non-uniform gases, 
Cambridge Univ. Press, Cambridge (1970). 

[8] J. Giacalone and J.R. Jokipii, ApJ 520 (1999) 204. 

[9] F. Casse, M. Lemoine and G. Pelletier, Phys. Rev. D 65 (2002) 023002. 

[10] J.R. Jokipii and E.N. Parker, ApJ 155 (1969) 777. 

[11] M.A. Forman, J.R. Jokipii, and A.J. Owens, ApJ 192 (1974) 535. 

[12] L. Spitzer Jr., Physics of Fully ionized gases, Interscience, New York (1962). 

[13] T. G. Northrop, The adiabatic motion of charged particles, Interscience, New 
York (1963). 

[14] B. Rossi and S. Olbert, Introduction to the Physics of Space, McGraw-Hill, New 
York (1970). 

[15] R.A. Burger, H. Moraal and CM. Webb, A&SS 116 (1985) 107. 

[16] G.K. Batchelor, The theory of homogeneous turbulence, Cambridge Univ. Press, 
Cambridge (1960). 

[17] J.W. Armstrong, J.M. Cordes, and B.J. Rickett, Nature 291 (1981) 561. 

[18] A. A. Ruzmaikin, A. M. Shukurov and D. D. Sokoloff, Magnetic fields of galax- 
ies, Kluwer Academic Press, Dordrecht (1988). 

[19] R.H. Kraichnan, Phys. Fluids 8 (1965) 1385. 

[20] A.M. Bykov and I.N. Toptygin, Ap&SS 138 (1987) 341. 

[21] R. Kubo, J. Phys. Soc. Japan, 12 (1957) 570. 

[22] M.A. Forman, Ap&SS 49 (1977) 83. 

[23] J.W. Bieber and W.H. Matthaeus, ApJ 485 (1997) 655. 

[24] L.J. Gleeson, Planet. Space Sci. 17 (1969) 31. 

16 



[25] P.A. Isenberg and J.R. Jokipii, ApJ 234 (1979) 746. 

[26] R. Balescu, H.-D. Wang, and J.H. Misguich, Phys. Plasmas 1 (1994) 3826. 

[27] V.S. Berezinsky et al., Astrophysics of Cosmic Rays, North- Holland, Amster- 
dam (1990). 

[28] D. Harari et al, JHEP 0203 (2002) 045. 

[29] G.G. Getmantsev, Sov. Astron. 6 (1963) 477. 

[30] I.H. Urch, Ap&SS 46 (1977) 389. 

[31] J. Kota and J.R. Jokipii, ApJ 531 (2000) 1067. 



17 



